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Abstract. Using ideal magnetohydrodynamics we examine an outflow from a disk surrounding a 
stellar-mass compact object. We demonstrate that the magnetic acceleration is efficient (J> 50% of 
the magnetic energy can be transformed into kinetic energy of 7 > 10 2 baryons) and also that the jet 
becomes collimated to very small opening angles. Observational implications, focusing on the case 
of an initially neutron-rich outflow, are discussed in Konigl's contribution. 



1. IDEAL MAGNETOHYDRODYNAMICS 

There is growing evidence in favor of magnetic driving in outflows associated with 
gamma-ray burst (GRB) sources [e.g., 5, hereafter VK03a; see also Konigl's contribu- 
tion in these Proceedings]. The dynamics of these outflows may be described to zeroth 
order by the ideal, axisymmetric, hydromagnetic equations, consisting of the Maxwell 
and momentum equations together with the conservations of baryonic mass and specific 
entropy. VK03a demonstrated that, under the assumptions of a quasi-steady poloidal 
magnetic field and of a highly relativistic poloidal velocity, these equations become ef- 
fectively time-independent and the motion can be described as a frozen pulse, generaliz- 
ing the so-called "frozen pulse" approximation already known in purely hydrodynamic 
models of GRB outflows [3]. Introducing the magnetic flux function A, the arclength 
along a poloidal streamline £, and the operator V 5 that acts while keeping s = ct — £ 
constant, the momentum equation can be written as (see VK03a for details) 

rpo( V.V s )( 4r V)=W±Miil-VP. (!) 

The large-scale electromagnetic field (E, B), the bulk flow speed (V), and the total (e ± + 
radiation) pressure can be written as functions of A and the rest baryon density po 1 

V S A x $ Q. AO. 2 ~ 4/3 

B= ^ +B , E= V ? A, V=- = — B + UJQQ, P = Qpl'\ (2) 

Faraday's law and the conservations of specific entropy and mass imply that the func- 
tions £1, Q, and Om are constants of motion, i.e., functions of A. By integrating equation 



(z,tD,<j)), and (r,0 ,<j>) denote cylindrical and spherical coordinates, whereas subscripts p and denote 
poloidal and azimuthal components, respectively. 



(1) along Y p and <j> one gets two additional constants of motion, 
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PLc 2 = ^yc 2 -^WB^ tt_A = ^ ymV ^_^ m ^ (3) 

describing the conservation of the ratio (total energy flux)/(mass flux) and of the total 
specific angular momentum, respectively. The remaining unknowns are the functions 
A(r,6) and Po(r, 6). The latter is the solution of Bernoulli's equation [which is obtained 
after substituting all quantities in the identity y 2 = 1 + y 2 V 2 /c 2 + y 2 V^ / c 2 using eqs. (2) 
and (3)], whereas the former controls the shape of the poloidal streamlines and is the 
solution of the highly nonlinear transfield component of the momentum equation (1). 

VK03a integrated the last two equations under the r self-similarity assumption A = 
r F /(0), which makes it possible (with Q. oc A~ l / F , Q<xA^ 2F ^ 3F , and constant F, Om> 
fl, and xa) to separate the (r, 6) coordinates. The resulting simplified set of ordinary 
differential equations can be easily integrated. Two types of boundary conditions at 
the base of the flow were considered: (1) The case of a strong poloidal magnetic field 
(Bp ^ BA, which corresponds to a trans- Alfvenic outflow (since the azimuthal field 
dominates asymptotically and the flow becomes super-fast). (2) The case of a strong 
azimuthal field (B p <B^), which corresponds to a super- Alfvenic flow. 



1.1. Trans-Alfvenic flows 

A representative solution is shown in Fig. 1. Looking at panel (a), which shows the 
acceleration, one can distinguish three different regimes: 

1) U5\ < (D < d6 is the fireball phase. The specific enthalpy £ decreases, resulting in 
increasing y oc gj y const, a characteristic of hydrodynamic acceleration), while 
the specific Poynting flux remains constant (the field is force free). The electromagnetic 
field only guides the flow, with the bulk of the collimation occurring in this regime. 

2) £J6 < £J < GJ8 is the magnetic acceleration regime. The fluid is cold (^ « 1), but y 
continues to increase (roughly as y oc jjj) due to the decreasing specific Poynting flux. 

3) GJ = U5% is the asymptotic cylindrical regime. The final Lorentz factor is of the order 
of the final specific Poynting flux, meaning that ~ 1/2 of the total energy (which was 
mostly electromagnetic initially) is transformed into baryonic kinetic energy (y^ w ju/2). 

The solution presented in Fig. 1 describes one shell, corresponding to a particular 
value of s. By specifying the s dependence in the initial conditions one can examine 
a multiple- shell outflow and the time dependence of the pulse. For example, the s 
dependence of the (total energy )/mass flux ratio }i(s)c 2 translates into different final 
Lorentz factors for distinct shells: f°°( s ) ~ In contrast with Michel's solution 

[1], in which the classical fast magnetosonic point is located at infinity and y^ w /i 1 / 3 , 
here this point is encountered at a finite height and most of the magnetic acceleration 
occurs further out, leading asymptotically to y^s) ~ /l(s) /2 ^> /ifs) 1 / 3 . Thus, not only 
is the magnetic acceleration highly efficient, but the stronger dependence of yoo(s) on the 
initial conditions [through fl(s)] can lead to a larger contrast in y«, between successive 
shells and hence to a higher efficiency of internal shocks [e.g., 2]. 




FIGURE 1. Trans-Alfvenic flow solution, (a) The Lorentz factor y, the ratio £, of the enthalpy to the 
rest energy, and the ratio of the Poynting flux to the rest-energy flux (top curve) are shown as functions 
of G7 along the innermost field line, (b) The meridional projections of the innermost and outermost field 
lines are shown on a logarithmic scale, along with a sketch of the central object/disk system. The field 
lines have a parabolic shape (z « G7 2 ) for G7 <, 10 9 cm and become asymptotically cylindrical. The vertical 
lines mark the positions of the various transition points along the innermost field line [see text and 4]. 



1.2. Super-Alfvenic flows 

A representative super-Alfvenic solution, corresponding to a base magnetic field 
(Btj, ~ 10 14 G, B p ~ \0~ 2 B$), is shown in Fig. 2. The super-Alfvenic solutions are 
distinguished from the trans-Alfvenic ones in two main respects: (1) During the ini- 
tial thermal-acceleration phase, some of the internal energy is transformed into electro- 
magnetic energy even as another part is used to increase V p . (2) During the subsequent 
magnetic-acceleration phase, the rate of increase of the Lorentz factor with z can be 
significantly lower than in the trans-Alfvenic case; the rate of increase of the jet radius 
with z is correspondingly higher. The overall magnetic-to-kinetic energy conversion ef- 
ficiency is higher. See [6] for further details and analytic scaling relations. 
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FIGURE 2. Same as Fig. 1, but for a super-Alfvenic solution. Here y °< GJ° 46 , z x GJ 1 '' 



Another potentially important aspect of super-Alfvenic outflows, namely, the possibil- 
ity that their initial composition is highly neutron-rich, could significantly alleviate the 
GRB baryon-loading problem. In [7] (see also Konigl's contribution) it is shown that, in 
contrast to the purely hydrodynamic case, the neutrons can decouple at a Lorentz factor 
that is over an order of magnitude smaller than for the protons. 
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